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ABSTRACT 

We use cosmological zoom-in simulations of galaxy formation in a Milky Way (MW)-sized halo started from 
identical initial conditions to investigate the evolution of galaxy sizes, baryon fractions, morphologies and an¬ 
gular momenta in runs with different parameters of the star formation-feedback cycle. Our fiducial model with 
a high local star formation efficiency, which results in efficient feedback, produces a realistic late-type galaxy 
that matches the evolution of basic properties of late-type galaxies: stellar mass, disk size, morphology dom¬ 
inated by a kinematically cold disk, stellar and gas surface density profiles, and specific angular momentum. 
We argue that feedback’s role in this success is twofold: (1) removal of low-angular momentum gas and (2) 
maintaining a low disk-to-halo mass fraction which suppresses disk instabilities that lead to angular momentum 
redistribution and a central concentration of baryons. However, our model with a low local star formation effi¬ 
ciency, but large energy input per supernova, chosen to produce a galaxy with a similar star formation history 
as our fiducial model, leads to a highly irregular galaxy with no kinematically cold component, overly extended 
stellar distribution and low angular momentum. This indicates that only when feedback is allowed to become 
vigorous via locally efficient star formation in dense cold gas, resulting galaxy sizes, gas/stellar surface density 
profiles and stellar disk angular momenta agree with observed z = 0 galaxies. 

Subject headings: cosmology: theory - galaxies: feedback - methods: numerical 


1. INTRODUCTION 

The Cold Dark Matter scenario, with a low mean matter 
density and a cosmological constant (ACDM), has proven to 
be broadly successful in explaining and predicting a variety 
of observations, such as the Cosmic Microwave Background 
temperature anisotropies (e.g Komatsu et al. 2011; Hinshaw 
et al. 2013; Planck Collaboration et al. 2014), the evolution of 
cluster abundance (Vikhlinin et al. 2009), and the large scale 
distribution of matter in the Universe (Conroy et al. 2006; 
Springel et al. 2006). Nevertheless, although the basic frame¬ 
work of galaxy formation within the CDM scenario (White & 
Rees 1978; Fall & Efstathiou 1980) is widely accepted, many 
aspects of the theory of galaxy formation are not yet fully un¬ 
derstood (for a recent review, see Silk & Mamon 2012). 

One of the most salient issues in galaxy formation mod¬ 
elling is understanding why galaxy formation is such an in¬ 
efficient process, i.e. what controls the low fraction of cos¬ 
mic baryons that end up as stars in the centers of dark matter 
haloes. A number of different methods, such as dark matter 
halo abundance matching (Tasitsiomi et al. 2004; Vale &. Os- 
triker 2004, 2006; Conroy et al. 2006; Conroy & Wechsler 
2009; Guo et al. 2010; Behroozi et al. 2013), satellite kine¬ 
matics (Klypin & Prada 2009; More et al. 2011), and weak 
lensing (Mandelbaum et al. 2006) (see Kravtsov et al. 2014, 
for a comprehensive discussion) point towards stellar to dark 
matter mass fractions of M*/Mh « 3-5% on average for L* 
galaxies and even smaller fractions for galaxies of smaller and 
larger mass. The fraction of baryons that ends up in observed 
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galaxies is thus well below the cosmological baryon fraction 
of Hb/Uni « 16% (Planck Collaboration et al. 2015). 

Such low baryon fractions are believed to be due to galactic 
winds driven by stellar feedback at the faint end of the stel¬ 
lar mass function (Dekel & Silk 1986; Efstathiou 2000) and 
by the active galactic nuclei (AGN) and the bright end (Silk 
& Rees 1998; Benson et al. 2003). Over the last two decades 
there has been an intense effort to incorporate and model these 
processes in cosmological simulations of galaxy formation 
(e.g., Katz 1992; Navarro & White 1993; Katz et al. 1996; 
Thacker & Couchman 2001; Stinson et al. 2006; Govemato 
et al. 2007; Scannapieco et al. 2008; Colin et al. 2010; Agertz 
et al. 2011; Avila-Reese et al. 2011; Guedes et al. 2011; Scan¬ 
napieco et al. 2012; Hopkins et al. 2011; Brook et al. 2012b; 
Agertz et al. 2013; Ceverino et al. 2014; Roskar et al. 2014; 
Booth et al. 2013; Aumer et al. 2013; Christensen et al. 2014b; 
Keller et al. 2014; Kimm et al. 2015; Murante et al. 2015). 
However, these processes have proven to be challenging to 
model and results are generally mixed. In particular, until 
recently simulations generally produced galaxies with larger 
than observed baryonic masses, and had difficulties in pro¬ 
ducing galaxies with realistic bulge-to-disk ratios (although 
see Governato et al. 2010; Guedes et al. 201 1). 

Recent work (e.g, Stinson et al. 2013a; Hopkins et al. 2014; 
Agertz & Kravtsov 2015) has demonstrated the importance of 
stellar feedback in not only predicting the z = 0 stellar mass- 
halo mass relation, but also its evolution, in particular during 
early stages of galaxy evolution (z > 2). Furthermore, in our 
previous study (Agertz & Kravtsov 2015) we have shown that 
the interplay between star formation and stellar feedback is 
complex, and multiple ways exists for simulations to repro¬ 
ducing global galaxy observables, such as stellar masses, flat 
rotation curves etc. However, internal characteristics of galax¬ 
ies, such as the Esfr - Ygas (Kennicutt-Schmidt) relation and 
the existence of thin galactic disks, are not reproduced by all 
models. 
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As modeling of star formation and stellar feedback in 
galaxy formation simulations improves, it is important to con¬ 
sider, and test against, a large set of galaxy properties. In par¬ 
ticular, it is widely recognized that in addition to basic proper¬ 
ties of galaxies, such as stellar mass and morphology, galaxy 
size and angular momentum play important roles in our un¬ 
derstanding of galaxy formation. 

In CDM models, angular momentum in galaxies is thought 
to originate during the initial phase of density perturbation 
growth, as collapsing peaks are tidally torqued by neigh¬ 
boring overdensities (Hoyle 1951; Peebles 1969; Doroshke- 
vich 1970; White 1984). Angular momentum of dark matter 
halos is often expressed as a dimensionless spin parameter 
A = y/\/2f?virPvir (as defined by Bullock et al. 2001), which 
characterizes the specific angular momentum in units of the 
angular momentum that would be required for rotational sup¬ 
port near the virial radius. Here f?vir and Vyir are the virial 
radius and virial velocity of the halo, and j the specific angu¬ 
lar momentum inside Ry^. 

In classical theories of disk galaxy formation (Fall & Efs- 
tathiou 1980; Ryden & Gunn 1987; Dalcanton et al. 1997; Mo 
et al. 1998), gas acquires the same specific angular momen¬ 
tum as the host dark matter halo, with only a fraction of it lost 
during gas condensation onto the halo center. The remaining 
angular momentum sets the disk size, and the fraction of the 
angular momentum that is ultimately retained is connected to 
galaxy morphology and the Hubble sequence (Fall 1983; Fall 
& Romanowsky 2013). 

While the spin parameter of the dark matter within the virial 
radius of halos is found to have log-normal distribution with 
a median A « 0.04 and rms variance of crinA ~ 0.55, regard¬ 
less of halo mass and cosmic time (e.g., Bullock et al. 2001; 
Vitvitska et al. 2002; Bett et al. 2010), recent studies have 
shown that accreted matter, dark matter as well as gas, upon 
entry of the virial radius, feature a higher spin than the halo’s 
average (e.g. Pichon et al. 2011; Kimm et al. 2011; Tillson 
et al. 2012; Danovich et al. 2012; Codis et al. 2012; Stew¬ 
art et al. 2013; Ubler et al. 2014; Danovich et al. 2015). In 
particular, Danovich et al. (2015) have recently used a set of 
zoom-in AMR simulations of galaxy formation at high red- 
shifts (z > 1 -2) to show that A of the cold gas when crossing 
the virial radius grows with time (see also Pichon et al. 201 1). 

Despite the large angular momentum content of accreted 
gas, numerical simulations of galaxy formation have long suf¬ 
fered from the “angular momentum catastrophe" (Navarro & 
White 1994; Navarro & Steinmetz 2000), whereby galaxies 
forming in the centers of dark matter halos are dominated by 
large bulges with a significantly lower specific angular mo¬ 
mentum content than their host haloes. In the last several 
years, galaxy formation simulations showed that this problem 
can be solved, or at least alleviated, by efficient stellar feed¬ 
back (Scannapieco et al. 2008; Zavala et al. 2008; Sales et al. 
2010). For example. Brook et al. (2011, see also Brook et al. 
2012a; Ubler et al. 2014; Christensen et al. 2014a) found that 
supernova feedback can selectively remove low angular mo¬ 
mentum gas via outflows, leading to disk formation. More 
recently, Genel et al. (2015) used the Illustris simulation suite 
to show that stellar and AGN feedback, tuned to reproduce 
the observed z = 0 stellar mass function, can produce a realis¬ 
tic distribution of specific angular momentum of galaxies. 

Despite the complexities of angular momentum evolution 
due to gas dynamics and effects of feedback outlined above, 
Kravtsov (2013) showed that observed galaxy sizes and ra¬ 


dial surface density profiles of baryons are strongly correlated 
with properties of their parent halo, and that sizes of both late 
type disks and early type spheroids appear to be set by the spe¬ 
cific angular momentum proportional to that acquired by the 
dark matter halo. Subsequent studies at higher redshifts show 
that this finding holds for star forming galaxies from z 8 to 
z R:! 0 (e.g., Shibuya et al. 2015, and references therein). This 
suggests that the angular momentum distribution in galaxy 
disks is quite similar to the predictions of the classical disk 
formation models (e.g.. Mo et al. 1998). 

There is a number of outstanding issues related to effects of 
stellar feedback in setting galaxy size and angular momentum. 
First, while feedback can prevent the loss of angular momen¬ 
tum content of galaxies, or remove the lowest angular momen¬ 
tum gas, hence promoting the formation of disk-dominated 
galaxies, it is not completely clear why the baryons remaining 
after feedback-driven outflows retain the angular momentum 
proportional to the average angular momentum of the halo. 
Second, as emphasized by Roskar et al. (2014), while strong 
feedback may help in explaining low batyon fractions, exces¬ 
sively violent energy and momentum input can lead to over¬ 
heating of gaseous disks, which leads to much thicker disk 
galaxies than are observed. At the very least, new feedback 
schemes, as they are developed and introduced, need to be 
tested against these basic empirical features of galaxy evolu¬ 
tion. 

In this work we use a suite of cosmological simulations of 
a Milky Way-mass halo started from the same initial condi¬ 
tions, but ran with different parameters for star formation and 
stellar feedback (presented in Agertz & Kravtsov 2015), to 
study the evolution of galaxy sizes, morphologies and angular 
momenta. Our goal is to 1) understand how different assump¬ 
tions about the efficiency of star formation and stellar feed¬ 
back affects the evolution of these fundamental properties of 
galaxies, and 2) compare these results to a variety of observed 
galaxy properties. 

We describe our cosmological simulations of galaxy forma¬ 
tion in § 2, present our results and compare them to observa¬ 
tional data in § 3, discuss our findings and compare to previ¬ 
ous studies on this subject in § 4, and summarize our results 
and conclusions in § 5. 

2. GALAXY FORMATION SIMULATIONS 
2.1. Star formation 

We carry out cosmological hydroH-A-body zoom-in simula¬ 
tions of Milky Way mass galaxies using the Adaptive Mesh 
Refinement (AMR) code RAMSES (Teyssier 2002). All simu¬ 
lations, and the star formation and feedback physics adopted, 
are presented in detail in Agertz et al. (2013) and Agertz & 
Kravtsov (2015). Briefly, we adopt a local star formation rate 
using the following equation: 

= (1) 

Iff 

where /h^ is the local mass fraction of molecular hydrogen 
(H 2 ), Pg is the gas density in a cell, % = y/3TTl32Gpg the lo¬ 
cal free-fall time of the gas, and Cff is the star formation effi¬ 
ciency per free-fall time. We adopt the model developed by 
Krumholz et al. (2008), Krumholz et al. (2009a), and Mc¬ 
Kee & Krumholz (2010), hereafter the KMT09 model, for 
the abundance of H 2 , based on radiative transfer calculations 
of idealized spherical giant atomic-molecular complexes sub¬ 
ject to a uniform and isotropic Lyman-Werner (LW) radiation 
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field. Relating star formation to the molecular gas is well mo¬ 
tivated empirically, as galactic star formation rate surface den¬ 
sities correlate well with the surface density of molecular gas, 
independent of metallicity, and poorly or not at all with the 
surface density of atomic gas measured on kpc scales (Bigiel 
et al. 2008; Krumholz et al. 2009b; Gnedin et al. 2009). 

We adopt star formation efficiencies per free fall time in the 
range of eg ^ 0.01 -0.1, motivated by observations of local 
GMCs (Lada et al. 2010; Murray 2011; Evans et al. 2014). 
As demonstrated by Agertz & Kravtsov (2015), even large 
local values of Cff (~ 0.1) can reproduce the low global star 
formation efficiency inferred from the Kennicutt-Schmidt re¬ 
lation (e.g. Bigiel et al. 2008) due to self-regulating effects of 
stellar feedback (see also Hopkins et al. 2011, 2014). 

2.2. Stellar Feedback 

We adopt the stellar feedback model described in Agertz 
et al. (2013). Briefly, each formed stellar particle is treated 
as a single-age stellar population with a Chabrier (2003) ini¬ 
tial mass function (IMF). We account for injection of energy, 
momentum, mass and heavy elements over time via super¬ 
nova type II (SNII) and type la (SNIa) explosions, stellar 
winds and radiation pressure (allowing for both single scatter¬ 
ing and multiple scattering events on dust) on the surrounding 
gas. Each mechanism depends on the stellar age, mass and 
gas/stellar metallicity, calibrated on the stellar evolution code 
STARBURST99 (Leitherer et al. 1999). Feedback is done con¬ 
tinuously at the appropriate times when each feedback process 
is known to operate, taking into account the lifetime of stars 
of different masses in a stellar population. To track the life¬ 
times of stars within the population we adopt the metallicity 
dependent age-mass relation of Raiteri et al. (1996). 

Momentum from stellar winds, radiation pressure, and SNe 
blastwaves is added to the 26 nearest cells surrounding a par¬ 
ent cell of the stellar particle. Thermal energy from shocked 
SNe and stellar wind ejecta is injected directly into the par¬ 
ent cell. We explore the concept of retaining some fraction 
of the thermal feedback energy in a separate gas energy vari¬ 
able over longer times than expected purely form the local gas 
cooling time scale. This approach was discussed by Agertz 
et al. (2013) and Teyssier et al. (2013) (for our choice of pa¬ 
rameters, see Agertz & Kravtsov 2015), and can be viewed 
as accounting for the effective pressure from a multiphase 
medium, where local unresolved pockets of hot gas exert work 
on the surrounding cold phase (see recent work on superbub¬ 
bles by Keller et al. 2014), or a placeholder for other sources 
of energy, such as turbulence and cosmic rays (Booth et al. 
2013). 

Heavy elements (metals) injected by supemovae and winds 
are advected as a passive scalar and are incorporated self- 
consistently in the cooling and heating routine. The code ac¬ 
counts for metallicity dependent cooling by using tabulated 
cooling functions of Sutherland & Dopita (1993) for gas tem¬ 
peratures 10^-10* ^K, and rates from Rosen & Bregman 
(1995) for cooling down to lower temperatures. Heating from 
the UV background (UVB) radiation is accounted for by us¬ 
ing the UVB model of Haardt & Madau (1996), assuming a 
reionization redshift of z = 8.5. We follow Agertz et al. (2009) 
and adopt an initial metallicity of Z = 1Q~^Zq in the high- 
resolution zoom-in region in order to account for enrichment 
from unresolved Pop III star formation (e.g. Wise et al. 2012). 

2.3. Simulation suite 


TABLE 1 

List of cosmological zoom-in simulations. 


Simulation'^ 

Description 


ALL_Efb_e001 

ALL_Efb_e001_5ESN 

ALL_Efb_e010 

All feedback, Cff = 1% 

All feedback, eff = 1%, ^snii = 5 x 10^^ erg 
All feedback, eff = 10% 


(a) All simulations adopt the KMT09 model (see main text), and second feed¬ 
back energy variable -Efb with /ft, = 0.5 and /^is = lOMyr. All simulations 
reach a minimum cell size of Ax = 75 pc. 


We adopt a WMAP5 (Komatsu et al. 2009) compatible 
ACDM cosmology with Ha = 0.73, Hm = 0.27, Hb = 0.045, 
tTg = 0.8 and Hq = 70kms~^Mpc“*. A pure dark matter sim¬ 
ulation was performed using a simulation cube of size Lbox = 
179Mpc. At z = 0, a halo of mass M 200 C ~ 9.7 x lO" Mq 
was selected for re-simulation at high resolution, and traced 
back to the initial redshift of z = 133. Here M 200 C is defined 
as the mass enclosed within a sphere with mean density 200 
times the critical density at the redshift of analysis. The cor¬ 
responding radius is r 2 ooc = 205 kpc. The mass within the ra¬ 
dius enclosing overdensity of 200 times the mean density is 
Mioom = 1.25 X 10^^ Mq and r 2 oom = 340kpc. When baryons 
are included in the simulations, the final total halo mass re¬ 
mains approximately the same. 

The selected halo does not experience any major merger 
after z 1.5, potentially favouring the formation of an ex¬ 
tended late-type galaxy. A nested hierarchy of initial con¬ 
ditions for the dark matter and baryons was generated using 
the GRAFIC-i-i-^ code, where we allow for the high resolu¬ 
tion particles to extend to three virial radii from the centre 
of the halo at z = 0. This avoids mixing of dark matter par¬ 
ticles with different masses in the inner parts of the domain. 
The dark matter particle mass in the high resolution region 
is ifiDM = 3.2 X 10^ Mq and the adaptive mesh is allowed to 
refine if a cell contains more than eight dark matter particles, 
and similar criterion is employed for the baryonic component. 
At the maximum level of refinement, the simulations reach a 
physical resolution of Ax « 75 pc. 

In this work we focus on three zoom-in simulations from 
the suite presented in Agertz & Kravtsov (2015). All 
models adopt the same feedback and star formation imple¬ 
mentation, but with diffferent parameters, and are picked 
to represent three different galaxy formation scenarios: 1) 
overcooling due to inefficient feedback (ALL_Efb_e001), 
2) efficient feedback via boosted supernova feedback 
(ALL_Efb_e001_5ESN), and 3) efficient feedback via effi¬ 
cient star formation (ALL_Efb_e010). 

In ALL_Efb_e001 the efficiency of star formation is 
eff = 1%. Agertz & Kravtsov (2015) demonstrated how this 
resulted in inefficient feedback incapable of driving global 
galactic outflows. This inability results leads to overproduc¬ 
tion of stars at high redshifts, in strong tension with semi- 
empirically derived star formation histories (e.g. Moster et al. 
2013; Behroozi et al. 2013). This run forms a very dense 
and rapidly rotating central concentration (see Section 3.5 and 
Figure 8 below), which makes the simulation computationally 
extensive. Given that the galaxy it produces was clearly unre¬ 
alistic, this simulation was run only to z = 1.5. 

For ALL_Efb_e001_5ESN the overcooling problem was 

^ http://graftc.sourceforge.net/ 
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resolved by boosting the released energy per supernova by a 
factor of five. This led to a good match to many considered 
observables at z > 1-2; semi-empirically derived star for¬ 
mation histories, the stellar mass-gas metallicity relation and 
evolution, the stellar mass-halo mass (M*-M 2 oo) relation and 
its evolution, as well as flat shapes of rotation curves. How¬ 
ever, the normalization of the Kennicutt-Schmidt relation was 
in tension with high redshift data, and the vigorous galactic 
outflows ended up preventing the formation of a cold galactic 
disk at any redshift, as we show below. 

ALL_Efb_e010 explores a scenario where vigorous 
galactic outflows are generated not by artificially boosting the 
available stellar feedback momentum or energy budget, but 
by adopting a higher local star formation efficiency per free- 
fall time, here eg = 10%, in agreement with observed values 
in massive GMCs (e.g. Lada et al. 2010; Murray 2011; Evans 
et al. 2014), which leads to more correlated stellar feedback 
events. This resulted in galactic properties in good match with 
all considered observables at z > 1, while producing galactic 
wind mass loading factors that decreases with increasing dark 
matter halo masses, leading to an epoch of disk formation at 
z < 1. In this work we continue the analysis down to z = 0. 

The simulations used in this study and their star formation 
and feedback parameters, are summarized in Table 1. 

3. RESULTS 
3.1. Morphology 

Eigure 1 shows mock composite g,r, i-fllter face-on and 
edge-on images of galaxies in the three simulations at their 
respective final redshifts, which account for dust attenuation. 
As we noted in the previous section, the simulation with low 
star formation efficiency (ALL_E f b_e0 01, bottom row) was 
not evolved to z = 0 due to the high computational expense, as 
the galaxy suffered from the classical overcooling problem 
(see § 3.2) and was stopped at z = 1.5. While the acquisition 
of galactic angular momentum, and diversity of galaxy mor¬ 
phologies, ultimately is linked to properties of cosmological 
flows and the halo merger history (see discussion in § 1), we 
here And that despite identical initial conditions, the resulting 
galactic morphologies are strikingly different due to relatively 
small changes in the parameters of small-scale star formation 
and feedback 

Overcooling in the ALL_E f b_e 001 run (bottom row) pro¬ 
duces a very massive galaxy and compact morphology with 
no clear signature of an extended stellar disk at any time. This 
morphology is in stark contrast with the well defined galac¬ 
tic disk formed in the fiducial model (ALL_Efb_e010, top 
row), which features both a thick and thin stellar disk com¬ 
ponent at z = 0. The vertical velocity dispersion of the young 
(f* < 3 Gyr) thin disk is cr*^ ~ lOkms”' (for a full analysis, 
see Agertz & Kravtsov 2015), and ~ 60-70kms“* for the 
old thick disk, compatible with what is observed in the Milky 
Way (e.g. Bovy et al. 2012). The model with low star for¬ 
mation efficiency, but boosted energy release per supernova 
(ALL_E f b_e 0 01_5ESN), fails to form a cold thin disk, with 
star formation occuring in an extended and highly turbulent 
gaseous system. The overall morphology is irregular, with 

* We note that we are only concerned with galaxies forming in ~ 
haloes, and not massive galaxies in the high-z Universe which likely are pro¬ 
genitors of today’s early type galaxy population (Zolotov et al. 2015). In 
massive disks at high z, properties of cold stream accretion, such as accretion 
alignment, likely matter more than for the galaxy types investigated here (e.g. 
Danovich et al. 2015) 



Weak feedback, z= 1.5 


15 kpc 


Fig. 1.— Mock g, r, i composite images of the three simulations at their fi¬ 
nal output, shown face-on in the left column and edge-on in the right column. 
The images are constructed assuming single stellar population SEDs from 
STARBURST99 (Leitherer et al. 1999) for each star particle, given its known 
age and metallicity, and then ray-tracing the line-of-sight flux. In the face-on 
images, the flux is attenuated using a MW-like reddening curve (Weingartner 
& Draine 2001) with a dust-to-gas ratio vs. gas metallicity relation compat¬ 
ible with Fisher et al. (2014). The galaxy in the run with low star formation 
efficiency (bottom row) features a very compact morphology at z = 1.5, in 
stark contrast to the fiducial model which features a well defined galactic 
disks at z = 0 (upper row). The model with boosted energy per supernovae 
(middle row) fails to form a thin disk, as star formation occurs in an extended 
and highly turbulent gaseous system. 


large gaseous clumps and irregular dark dust lanes. 

3.2. Baryonfractions 

Eigure 2 shows the stellar mass fraction of the central galax¬ 
ies in the three runs, expressed in units of the cosmic baryon 
fraction, as a function of dark matter halo virial mass for the 
three models. We compare the simulation data to results from 
semi-empirical constraints based on the abundance matching 
approach (Behroozi et al. 2013; Kravtsov et al. 2014), as well 
as inferences from weak lensing studies of late and early type 
systems (Mandelbaum et al. 2006; Hudson et al. 2015) and 
estimates using satellite kinematics (More et al. 201 1). 

The fiducial and the boosted feedback runs conform to the 
z = 0 stellar mass - halo mass relation, and its inferred evolu¬ 
tion. As already discussed in Agertz & Kravtsov (2015), these 
models achieve this in very different ways leading to dramatic 
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Fig. 2.— The evolution of the stellar mass fraction in units of the universal 
baryon fraction as a function of halo mass. The shaded regions show, from 
light to dark grey, the z = 3,2,1 data from Behroozi et al. (2013) and z = 
0 from Kravtsov et al. (2014), where the gray bands show the region ±2cr 
around the mean relation. The points show observational inferences of stellar 
fractions for galaxies at z 0 using weak lensing (circles) for late type (blue) 
and early type (red) systems by Mandelbaum et al. (2006) and Hudson et al. 
(2015). Magenta triangles show estimates using satellite kinematics by More 
et al. (2011). 

differences in the resulting galaxy morphologies (Figure 1). 
The run with low local star formation efficiency is inconsis¬ 
tent with high redshift data, and is marginally consistent with 
the expected scatter of stellar fractions at z = 0 , although note 
that the hnal data point is at z = 1.5. 

3.3. Size evolution in simulations and comparison with 
observations 

Figure 3 shows evolution of the stellar mass-half mass ra¬ 
dius (Mi,-r\i'^ relation for galaxies in our three runs against 
the observed relation at z = 0 (see Kravtsov 2013, for details 
about observational sample and the method used to estimate 
ri/ 2 ). The low efficiency star formation run evolves off the 
mean relation at early times, and at the last available data point 
(z= 1.5) the galaxy is extremely compact with ri /2 ^ 0.5 kpc. 
We note that a signihcant fraction of z 1-3 galaxies are 
in fact compact ellipticals (van Dokkum et al. 2008; Newman 
et al. 2010; van Dokkum et al. 2014) with r^ji < Ikpc (Car- 
ollo et al. 2013). However, such galaxies are believed to be 
progenitors of today’s elliptical galaxies that evolve at low z 
primarily due to mass accretion via minor mergers (e.g., Hop¬ 
kins et al. 2009). The compact ellipticals are mostly quenched 
at low z, as characterized by their specihc star formation rates. 
With a SFR^ 2QMq yr“* at z = 1.5, the overcooling model is 
far from being quenched, and while it could be identihed as 
a “blue nugget" (progenitors of compact ellipticals, see e.g. 
Barro et al. 2013; Dekel & Burkert 2014), and join the pop¬ 
ulation of early type galaxies at z = 0 by quenching and late 
time expansion (Whitaker et al. 2012), its already signihcant 
stellar mass fraction makes it a rather unlikely candidate for a 
realistic galaxy. 

In contrast, the boosted feedback run 
ALL_Efb_e001_5ESN, in which the resulting galaxy 
lacks a cold disk component, is very extended, with 
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Fig. 3.— Evolution of the stellar 3D mass-half mass radius relation for 
galaxies in our three runs (lines) against the observed relation at z = 0 (points) 
for galaxies of different morphological types. The simulations are shown by 
the solid blue (fiducial), magenta (ALL_Efb_e001_5ESN ), and dashed 
light blue line (low star formation efficiency, weak feedback) lines. Ob¬ 
served galaxies shown by points of the following type: the red pentagons 
and hexagons show a sample of elliptical and dwarf elliptical galaxies from 
the compilation of Misgeld & Hilker (2011); blue circles are the late type 
galaxies from the samples of Leroy et al. (2008) and Zhang et al. (2012) with 
half-mass radii estimated from the deprojected surface density profiles, while 
the star symbol shows the Milky Way; the red cartwheel points show the 
Local Group dwarf spheroidal galaxies from the compilation of Misgeld & 
Hilker (2011). The light blue dashed line and dot-dashed orange line show 
the average relations derived by Bemardi et al. (2014) for late and early-type 
galaxies, respectively. Dark red shaded band shows 2cr scatter around the 
mean relation calculated for all galaxies in the Bemardi et al. (2014) sample. 
Individual galaxies from the SDSS, presented in Szomoru et al. (2013), are 
shown by dots. 


ri /2 > lOkpc at z = 0. This is ^2-3 times greater than the 
observed average ri/ 2 - Interestingly, extended spheroidal 
systems with exponential surface brightness profiles do exist 
in the nearby universe (van Dokkum et al. 2015). However, 
such galaxies are rather rare and cannot be considered to be 
the typical outcome of galaxy formation in a « 10 '^ Mq halo. 

Figure 3 also shows that the galaxy forming in the fiducial 
run matches the z = 0 observations at all times. Prior to the last 
major merger at z ~ 1 . 8 , when M* « 10 '® Mq, ri/2 is roughly 
constant at ~ 1-2 kpc. After this point the galaxy enters an 
epoch of disk formation and grows rapidly in size, with rj /2 ~ 
3-4 kpc at z = 0, close to the estimated value for the Milky 
Way (Bovy & Rix 2013). More precisely, starting fromM* « 
1 O*®M 0 at z 1 . 8 , our fiducial model grows by a factor of ~ 
1.75 in ri /2 while growing by a factor of 6 in mass until z = 0, 
indicating rj /2 cic M® 3. This is in excellent agreement with the 
the results of van Dokkum et al. (2013) and Patel et al. (2013), 
who studied the structural evolution of MW-mass progenitors 
back to z ~ 2.5 using the 3D-HST and CANDELS Treasury 
surveys. Here the half-light radii of star forming galaxies in 
their sample were found to scales as ri /2 c>c ^0-29±o.08 

These results illustrate that while the galactic environment, 
via mergers, alignment of accretion etc., must play a role 
for the resulting galaxy size and morphology (especially in 
explaining the origin of early type galaxies, e.g., Hopkins 
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Fig. 4.— Left: 3D half-mass radius- R 200 radius relation of galaxies, where S 200 is the radius enclosing density contrast of 200 relative to the critical density. 
Right: the same relation but with i ?200 of observed galaxies corrected approximately for pseudo-evolution (see text for details). The simulations are shown by the 
solid blue (fiducial), magenta (ALL_Efb_e001_5ESN ), and dashed light blue line (low star formation efficiency, weak feedback) lines, which are identical in 
both panels. The /? 20 o(z) radius is measured for the simulated galaxies at the redshift of analysis. As shown by Kravtsov (2013), the observed galaxies at z Ri 0 
on average follow a nearly linear relation, r [^2 ~ O.OI 5 R 200 spanning more than eight orders of magnitude in stellar mass. The red pentagons and hexagons show 
a sample of elliptical and dwarf elliptical galaxies from the compilation of Misgeld & Hilker (201 1); blue circles are the late type galaxies from the samples of 
Leroy et al. (2008) and Zhang et al. (2012) with half-mass radii calculated directly from the deprojected surface density profiles, while the star symbol shows 
the Milky Way; the red cartwheel points show the Local Group dwarf spheroidal galaxies from the compilation of Misgeld & Hilker (2011). The light blue 
dashed line and dot-dashed orange line show the average relations derived for late and early-type galaxies, respectively, from the average R 112 — M* relations of 
Bernardi et al. (2014). Dark red shaded band shows 2a scatter around the mean relation calculated for all galaxies in the Bemardi et al. (2014) sample. The 
orange dot-dashed line with error bars shows the mean relation and 2a scatter for massive SDSS galaxies presented in Szomoru et al. (2013); individual galaxies 
from this sample are shown by blue (Sersic index n < 2.5) and red {n > 2.5) dots. The gray dashed line shows linear relation riy 2 = 0.015i?200 and dotted lines 
are linear relations offset by 0.5 dex, which approximately corresponds to the scatter in galaxy sizes from distribution of halo spin parameter A under assumption 
that ri /2 oc Ai? 200 - 
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et al. 2010), internal small-scale feedback processes, and our 
limited theoretical understanding of how these should be ac¬ 
counted for over cosmic time, impacts, and possibly even dic¬ 
tates, the Hubble sequence of late type galaxies. We discuss 
this further in § 4. 

The left panel of Figure 4 shows the evolution of galaxies in 
the three runs in the plane of half-mass radius, ri/ 2 , vs. the ra¬ 
dius, Rioo, enclosing the density contrast of 200 relative to the 
critical density of the universe. Simulation results are shown 
against a collection of z « 0 observed galaxies (see caption 
for the data sources), for which /?200 was estimated from the 
halo mass obtained by abundance matching. As shown by 
Kravtsov (2013), this relation is on average close to linear, 
ri /2 « 0.015/? 200 , which is shown by the long gray dashed 
line along with the 2a scatter, as expected if galaxy sizes were 
set by the specific angular momentum of baryons shared with 
the dark matter halo (i.e., ri /2 ex Af? 2 oo)- A similar relation has 
since been observed at a wide range of redshifts (e.g., Shibuya 
et al. 2015). 

The figure shows that in the weak feedback simulation the 
galaxy maintains a compact distribution of stars with half¬ 
mass radius staying constant at ri /2 ~ 0.3-0.5 kpc. The 
fiducial and ALL_Efb_e001_5ESN simulations, however, 
evolve approximately along the linear relation, albeit with 
significant fluctuations, especially during early stages evolu¬ 
tion (smaller sizes). As discussed above, the galaxy in the 


ALL_Efb_e001_5ESN run has a size that is too large for 
a typical halo expected to host its galaxy at z = 0. Such galax¬ 
ies are quite rare in the observed samples. The fiducial run is 
close to the z = 0 data inferred for observed galaxies during 
the late stages of evolution. At high z, the ratio r 1 / 2/^200 is 
about twice larger than at low z, a behavior that is consistent 
with the average behavior of observed galaxies reported by 
Shibuya et al. (2015). Such evolution may be partly due to 
the physical evolution of the disk size, but could also arise be¬ 
cause the disk size is set during the fast mass growth regime of 
halo evolution after which a large fraction of the M 200 growth 
is not a real physical growth but an increase of mass due to 
pseudo-evolution (Diemer et al. 2013), and is hence not ac¬ 
companied by the proportional growth of stellar ri/ 2 - 
To illustrate the latter point, the right panel of Figure 4 
shows a somewhat different representation of the evolution. 
Here, the curves representing tracks of simulated galaxies are 
unchanged from the left panel. The observed z = 0 sample, 
however, is corrected for possible pseudo-evolution (PE) of 
the halo mass. Many indications exist that the inner regions 
of halos change little from z 1 - 2 until the present (Cuesta 
et al. 2008; Diemer et al. 2013; More et al. 2015). Thus, the 
evolution of M 200 and R 200 at these redshifts may be sim¬ 
ply due to the changing critical density of the universe, to 
which they are tied by definition, and not by the real phys¬ 
ical growth. In the extreme scenario, one can imagine that 








Feedback and galaxy structure 


7 



R (kpc) 



R (kpc) 


Fig. 5.— Comparison of the surface density profiles of stars (left) and cold gas (right) in the simulated galaxies in the fiducial and ALL_Efb_e001_5ESN 
simulations along with the profiles of late type galaxies from the THINGS sample (Leroy et al. 2008) in the stellar mass range M* ps 3 — 10 x 10*^ Mq and 
the Milky Way, using a combination of the thin and thick stellar disks with parameters given in Table 2 of McMillan (201 1) for S* and the sum of HI and H 2 
profiles from Dame (1993) for Sgas- In observed galaxies the gas surface density includes atomic and molecular gas corrected for Helium mass fraction, while in 
simulations the gas mass includes all of the cold gas of T < lO'^ K. 


halo grows only by pseudo-evolution since its formation time, 
which we can operationally define as the time when the halo 
had a concentration of Cform = 3 (Wechsler et al. 2002; Zhao 
et al. 2009). Thus, in the right panel we correct the virial radii 
as = Cform7?200/c(M200,z = 0), whcrc c(M200,z = 0) is the 
average concentration of halos of mass M 200 at z = 0, which we 
calculated using an accurate concentration model by Diemer 
&. Kravtsov (2015) calibrated on a large suite of simulations. 

If galaxy halos do exhibit signihcant pseudo-evolution, the 
disk sizes were set at higher z and thus they should be com¬ 
pared to the PE-corrected /?200 values. In other words, in such 
case it should be more appropriate to compare higher-z por¬ 
tion of the galaxy tracks with the data with the PE-corrected 
virial radii. The right panel of Eigure 4 indeed shows that the 
early stages of galaxy evolution, and hence lower progenitor 
stellar masses and halo virial radii, track the PE-corrected data 
points. At late stages they evolve close to the z « 0 data. This 
evolution qualitatively and quantitatively tracks the evolution 
derived for observed galaxies across a wide range of redshifts 
by Shibuya et al. (2015), who find that star forming galaxies 
have ri/ 2/^200 ~ 0.015-0.02 ?A.z<2 increasing somewhat to 
'" 1 / 2/^200 ~ 0.03 - 0.04 at larger redshifts. 

Overall, these results indicate that in the fiducial and 
ALL_Efb_e001_5ESN simulations, the angular momen¬ 
tum of the disk is approximately linearly proportional to the 
specific angular momentum of the halo and galaxies thus 
evolve along the ri /2 oc /? 2 oo track. Similar results was re¬ 
cently reported by Pedrosa & Tissera (2015), who found 
that galaxies in a large simulated sample followed a linear 
^200 “^ 1/2 relation in simulation with efficient feedback, but 
disk and spheroidal galaxies followed relations with different 
constants of proportionality. 


3.4. Surface density profiles 

Eigure 5 shows the stellar and neutral gas surface den¬ 
sity profiles for the simulated galaxies in the fiducial and 
ALL_Efb_e001_5ESN simulations - the two simulations 
that were evolved to z = 0 - along with the profiles of late type 
galaxies from the THINGS sample (Leroy et al. 2008) in the 
stellar mass range M* « 3 -10 x 10'° Mq and the Milky Way 
(McMillan 2011; Dame 1993). In the observed galaxies the 
gas surface density includes atomic and molecular gas cor¬ 
rected for Helium mass fractions, while in the simulations the 
gas mass includes all of the cold gas of T < 10"^ K. 

The figure shows that the fiducial simulation reproduces 
stellar surface density profile of late type galaxies of this mass 
range, albeit with slightly larger bulge component than typi¬ 
cal in the observed sample. Thus, this simulation does not 
just result in a stellar distribution with the correct half-mass 
radius, but also naturally reproduces a radial prohle compati¬ 
ble with observed stellar distributions. In contrast, the surface 
density prohle in the ALL_Efb_eO 01_5ESN simulation is 
very different from all of the observed THINGS galaxies in 
this stellar mass range; the E*(7?) prohle is close to exponen¬ 
tial, but as discussed above the scale length is considerably 
larger than most disk galaxies. Thus, it is a counterpart to 
the rare low-surface brightness galaxies of which some fea¬ 
ture very extended exponential disk (e.g., McGaugh & Bothun 
1994). Nevertheless, the number density of such galaxies, 
is roughly two orders of magnitude smaller than the num¬ 
ber density of typical L* galaxies we are investigating. Such 
galaxies are thus unlikely to form in a single random realiza¬ 
tion of a MW-sized halo. 

The gas surface density prohle of the hducial simulation 
does not match the prohle of the Milky Way. However, its 
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Fig. 6.— Evolution of the specific angular momentum of galaxies in the three studied simulations and comparisons with specific angular momenta of galaxies 
of different moiphological type from Fall & Romanowsky (2013). The panels show stars (left), cold T < 10“* K gas (middle), and the total specific angular 
momentum of stars and cold gas (right). 


shape matches the overall shape of the gas surface density 
profiles in this mass range. The normalization is somewhat 
higher than most THINGS galaxies which indicates that the 
simulated galaxy may be more gas rich than typically ob¬ 
served, and we return to the implications of this in § 4. 

The Sgas(r) profile in the ALL_Efb_e001_5ESN simu¬ 
lation has an even higher normalization. Overall, this galaxy 
has lower star formation rate than the fiducial simulation at 
low redshifts, 0.5 < z < 2 (Agertz & Kravtsov 2015) and it 
retains more gas. The shape of the gas surface density pro¬ 
file is somewhat steeper than the prohles of observed galax¬ 
ies. Profiles in both simulations exhibit more fluctuations than 
the due to the presence of spiral arms and other in¬ 

homogeneities in the gas distribution. Although the distri¬ 
bution of young stars exhibits similar fluctuations, the stellar 
mass distribution is considerably smoother which is reflected 
in the smoother profile. Overall, radial fluctuations of 

I]gas(I?) in simulated galaxies are comparable to those in the 
data, which indicates a similar level of structure. 

3.5. Evolution and scaling of angular momentum 

Figure 6 shows the evolution of stellar, cold gas {T < 10^ K) 
and joint specific angular momentum^, of the central galaxies 
in the three simulations as a function of stellar mass. The 
models are compared to specific angular momenta of galaxies 
of different morphological types from Fall & Romanowsky 
(2013), where the authors estimate the disk-only contribution 
via i = IVmtRd, where is the scale radius of the exponential 
disc and V^ot is the asymptotic rotation velocity observed at 
large radii. 

The M 200 “ M* and M* - rj /2 relations in Figure 2 and 3 
indicated a monotonic relation between feedback strength and 
final galaxy size; as the strength of feedback increases, ri /2 
increases and M* is decreased at any given value of M 2 oo- We 
find that the same simple behavior is not recovered for the 
stellar mass - specific angular momentum relation (M*-J), as 
clearly visible for all ISM components in Figure 6 . 

The galaxy in the fiducial run has a specific stellar angu¬ 
lar momentum 7 * < 100 kms“'kpc during early stages when 
M* < IO'^Mq (z > 1.5). Angular momentum before that 

^ The specific angulai' momentum is computed as j = 'ti/ r,- ’^i 

where summation is performed over all i cells of cold gas or all star particles, 
and r and v are position and velocity relative to the galaxy’s center of mass. 


epoch fluctuates due to repeated mergers with the main pro¬ 
genitor leading to bursts of star formation and vigorous out¬ 
flows (see Agertz & Kravtsov 2015, for a thorough discus¬ 
sion). The last major merger occurs around z ^ 1.8, after 
which an epoch of disk formation begins, correlating with an 
increase in 7 * by almost an order of magnitude. The corre¬ 
sponding increase in ri /2 is only a factor of ^ 2 as discussed 
above. The final specific angular momentum of the galaxy in 
the fiducial run is compatible with spiral galaxies of roughly 
Sa-Sab type. 

In the run with inefficient feedback, 7 * is always < 
100 kms“*kpc, in agreement with early results on the “an¬ 
gular momentum problem" (e.g. Navarro & Steinmetz 2000). 
This does not come as a surprise, as the specific angular mo¬ 
mentum content of accreting gas is expected to be smaller at 
higher redshifts, for a given galaxy progenitor (e.g. Pichon 
et al. 2011; Danovich et al. 2015), and this galaxy formation 
model overpredicts the fraction of stellar mass forming out of 
this gas at early times instead of ejecting it via galactic winds 
(see also Ubler et al. 2014). In addition, the fraction of mass 
in the central galaxy is above the stability threshold for disks 
against an m = 2 instability and bar formation. Such instabil¬ 
ities channel the gas towards the center, where it forms stars 
leading to a lower ri/ 2 - We discuss this issue in more detail in 
Section 4.1 below. 

Rather surprisingly, the boosted feedback model fails to 
produce a galaxy rich in angular momentum, with 7 * ~ 100 - 
500kms“’kpc, with large scatter, despite the large physical 
extent as characterized by ri/ 2 - In contrast to the effect of 
feedback on the M*-ri /2 evolution, the interaction between 
feedback and specific stellar angular momentum is more com¬ 
plex. In the case of the ALL_Efb_e001_5ESN run, strong 
feedback results in a more spheroidal distribution of stars with 
angular momenta of individual orbits partially canceling each 
other. 

These results imply that strong feedback is not a panacea 
against the angular momentum catastrophe, as different mod¬ 
els of the star formation - feedback cycle can give almost 
identical results in terms of global quantities (M*, SFR, Vdrc), 
but differ significantly in other, more detailed, observables 
(Agertz & Kravtsov 2015). In §4.2 we discuss feedback im¬ 
plementations further, and compare our results to recent stud¬ 
ies on the angular momentum content of simulated galaxies. 
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Fig. 7.— The evolution of the halo baryon mass fraction (stars and cold 
gas) in the central galaxy in the three simulations (the colors and types ai‘e 
the same as in Figure 6). The shaded region indicates a transition where the 
disk is expected to be susceptible to bar formation, leading to mass transport 
to the center (see text for detailed discussion). 


We find that the angular momentum in the cold gas com¬ 
ponent (y'gas) is always greater than the stellar counterpart. 
This is not surprising as stars form in the densest inner re¬ 
gions of the gaseous disk. Without stellar feedback, y'gas 
reaches ~ 1000kms“'kpc at z ~ 1.5, 10 times higher than 
the corresponding y*. In the boosted feedback case, y'gas ~ 
500kms“'kpc after the last major merger (i.e. for M* > 
IO'^’M©), with significant scatter. The fact that jgas here is 
smaller than that in the fiducial model, by a factor of a few, 
may stem from how stellar feedback driven outflows, in this 
particular model, not preferentially ejects angular momentum 
poor gas (as demonstrated by Brook et al. 2011; Ubler et al. 
2014), but also high angular momentum gas, hence destroy¬ 
ing the cold ISM of the galaxy; by artificially boosting the 
available energy per SN, even gas associated with regions of 
low Esfr, which otherwise would not necessarily be as af¬ 
fected by feedback, can be ejected from the galaxy. Overall, 
in the runs with feedback jgas evolves approximately along the 
j-Mif relation of observed galaxies. 

In the fiducial model, jgas 2000km s“'kpc, a value com¬ 
patible with observed y* in late type spirals of Sc-Sd type. We 
note that the precise value for the simulated stellar angular 
momentum is sensitive to numerically resolving star forma¬ 
tion at large radii. Furthermore, the cold gas extends much 
farther than the optical radius (see Figure 5), in agreement 
with observations (e.g., Kravtsov 2013). 

4. DISCUSSION 

The results presented in the previous section clearly show 
that the morphology of a galaxy forming in a given dark mat¬ 
ter halo sensitively depends on the parameters of star forma¬ 
tion and feedback. This conclusion is in agreement with a 
number of other recent studies, as we discuss below in Sec¬ 
tion 4.2. 

Although the role of feedback in setting galaxy morphology 
and basic properties, such as stellar mass and size, is widely 


acknowledged, the exact mechanisms by which feedback af¬ 
fects galaxy properties are not yet fully understood. Our re¬ 
sults can be used to clarify some aspects of this question, as 
we discuss next. 


4.1. Disk instabilities, bulge formation and stellar feedback 


It is widely accepted that stellar and AGN feedback play 
key role in removing low angular momentum gas that is ex¬ 
pected in the central regions of galaxies, even if there is no 
re-distribution of angular momentum during galaxy forma¬ 
tion (e.g., Bullock et al. 2001; van den Bosch 2001; van den 
Bosch et al. 2001). Flowever, the role of feedback is likely 
to be more multi-faceted. Even if cosmologically accreted 
gas has initially no low angular momentum component, re¬ 
distribution of angular momentum due to disk instabilities can 
quickly create such component (e.g., Kormendy & Kennicutt 
2004; Bournaud et al. 2007; Elme^een et al. 2008; Wein¬ 
berg & Katz 2007). If star formation proceeds in this low 
angular momentum gas, this results in a more massive central 
bulge and possibly a spheroid-dominated galaxy. If, on the 
other hand, stellar feedback removes most of the gas before 
it is turned into stars, the remaining disk may have a much 
smaller stellar mass much than in observed galaxies. The role 
of feedback should thus be to remove low angular momentum 
gas efficiently but still reproduce observed disk masses. As 
we argue below, such self-regulating behavior can probably 
be understood by considering overall disk stability. 

Disk stability is affected by a number of factors (e.g. Jog 
& Solomon 1984; Romeo 1992; Christodoulou et al. 1995; 
Rafikov 2001; Romeo et al. 2010; Agertz et al. 2015), but a 
key factor controlling the overall global stability is the the 
fraction of potential contributed by the disk relative to the 
spheroidal component (Efstathiou et al. 1982). This depen¬ 
dence was realized by ©striker & Peebles (1973), who used it 
to argue that existence of dynamically cold galactic disks re¬ 
quires the existence of massive spheroidal halos around them. 

Efstathiou et al. (1982) used A-body simulations to show 
that for exponential discs embedded in spherical haloes the 
onset of the bar instability can be characterized by the crite¬ 
rion* 


Et, 


-\/GMdisk/Rd 


< 1 - 1 , 


( 2 ) 


where Vmax is the maximum rotational velocity, Mdisk is the 
mass of the disk, and the disk’s exponential scale length. 
Due to its simplicity, this is a widely used criterion of insta¬ 
bility in semi-analytical models of galaxy formation (e.g. Mo 
et al. 1998; De Lucia et al. 2011; Guo et al. 201 1). 

Mo et al. (1998) used a detailed disk formation model to 
show that the criterion in Equation 2 implies disk instabil¬ 
ity when the disk baryon mass relative to the total halo mass, 
md = Md/M 2 oo, is greater than the disk spin parameter Ad: 
'Wd ^ Ad. If the disk has a specific angular momentum sim¬ 
ilar to that of the dark matter halo, its spin parameter is ex¬ 
pected to have a log-normal distribution with the median of 
A « 0.04 with a fractional scatter of crinA ~ 0.55, regardless 
of halo mass and cosmic time (e.g., Bullock et al. 2001 ; Vitvit- 
ska et al. 2002; Bett et al. 2010). This result implies that 
when galaxy growth increases the disk mass beyond this limit, 
strong disk instabilities will set in. 


^ It is worth keeping in mind that other factors such as the vertical stellar 
velocity dispersion and multiple disk components will affect this criterion, as 
pointed out by Athanassoula (2008). Thus, this criterion is approximate. 
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Fig. 8.— Circular velocities (Vc = \/GMjr) for the three galaxy models at j = 3 (left), 1.5 (middle) and 0 (right). The thick lines show Vc derived for all mass, 
whereas the thin lines account only for the stars and cold gas situated in the disk. 


The subsequent evolution will depend on the efficiency of 
feedback; if feedback removes the low angular momentum 
gas resulting from momentum re-distribution in the unstable 
disk, and lowers the disk mass below the critical limit, the 
galaxy may maintain a disk mass fraction close to the critical 
instability threshold. If feedback is inefficient, the galaxy will 
form a massive central spheroidal component. 

These different scenarios are indeed realized in our runs 
with different star formation and feedback parameters. In Fig¬ 
ure 7 we plot the disk baryon fractions (stellar -i- cold gas) for 
the three simulated galaxies as a function of redshift, with 
the grey band roughly indicating a possible transition into 
a disk unstable regime, and hence bar formation, assuming 
A = j/V 2 R 200 V 200 = 0.044, which we measure for our halo at 
z = 0. Note that the mass fraction includes all of the mass in 
the galaxy, while stability criterion should include only disk 
mass. We do not attempt to separate disk and spheroidal com¬ 
ponents and simply examine the total mass of the galaxy. The 
rationale is that if the entire mass of the galaxy was in a cold 
disk, the figure would shows us the regime in which the disk 
is expected to be unstable. 

The run with low efficiency of star formation does not 
have sufficiently strong feedback and produces a galaxy with 
^ 80% of the cosmic baryon fraction in the disk at all times. 
The feedback in the ALL_Efb_e001_5ESN run, however, 
does drive strong outflows at z > 1 and maintains disk masses 
right around the critical instability threshold from z 5 to 
z = 0. The feedback in the fiducial run is even stronger at 
high z, which keeps the disk mass is well below the insta¬ 
bility threshold until the last major merger at z ~ 1.8. At 
lower z, outflows become less prominent (Agertz & Kravtsov 
2015; Muratov et al. 2015) and the disk mass fraction reaches 
ma = Mdi.sk/M 200 ~ 0.07 (or « 0.4-0.5 in units of flb/f^m) at 
lower z. This could indicate that the fiducial model is suscep¬ 
tible to m = 2 disk instabilities at late times, which are respon¬ 
sible for growth of its bulge. 

Indeed, this is consistent with evolution of mass distribu¬ 
tion traced by the circular velocity profile in Figure 8, which 
shows the circular velocity for all three runs at different z for 
only the baryons in the disk and all matter separately. At 
z = 3, ALL_Efb_e001_5ESN and the fiducial runs both 
have disk masses below the critical instability threshold and 
correspondingly have circular velocity profiles without a cen¬ 
tral “spike”, which is very pronounced in the profile for 


the weak feedback run. At z = 0, on the other hand, the 
Vc(r) profile for the fiducial model features a modest peak in 
Vc(r) at r < 1 kpc consistent with late evolution of this disk 
above the critical threshold, whereas the circular velocity in 
the ALL_Efb_e001_5ESN run, in which the galaxy mass 
never gets significantly larger than the threshold, maintains 
a slowly rising Vc(r) profile without a central concentration. 
This figure thus indicates that development of central concen¬ 
trations in galaxies is related to the development of disk insta¬ 
bilities when the central disk mass exceeds the approximate 
threshold of Equation 2 above. 

This sheds light on the long-standing issue of “spiky” cir¬ 
cular velocity profiles in simulations (e.g., Mayer et al. 2008, 
and references therein): our results indicate that prevalent 
central spikes in Vc(r) profiles, even at the highest z, are due to 
the inability of a given feedback prescription of keeping the 
disk mass around or below the instability threshold. 

At the same time, our results show that preventing forma¬ 
tion of the central concentration does not necessarily result in 
formation of a disk-dominated galaxy. The galaxy that forms 
in our ALL_Efb_e001_5ESN run has maintained low disk 
mass and has exponential surface density profile, but has an 
irregular morphology and a disk size that is larger than ex¬ 
pected for the galaxy of such stellar mass. This is because 
stellar feedback in this run continues to strongly “stir” the gas 
in the disk and disrupt galaxy ISM with outflows and foun¬ 
tains down to the smallest redshifts. The ISM in this run 
is thus perpetually in a perturbed state, which prevents for¬ 
mation of a well-defined cold disk component (Roskar et al. 
2014). 

The above analysis illustrates the dual effect of stellar feed¬ 
back: it both removes the low-angular momentum gas and 
regulates disk stability. In late type L ^ L* galaxies, i.e. galax¬ 
ies with low bulge-to-disk-ratios or lack of bulges altogether, 
are predominately low mass galaxies (M* < 1-5 X 10'°Mo 
Nair & Abraham 2010), feedback effects should delicately 
balance between driving ISM turbulence and outflows and al¬ 
lowing thin disks to form. Such galaxies likely form in haloes 
of masses ^ 10'' - IO'^Mq (Behroozi et al. 2013; Kravtsov 
et al. 2014). 

van Dokkum et al. (2013) studied the structural evolution 
of MW-mass progenitors back to z ~ 2.5 by matching cumu¬ 
lative co-moving number densities in the 3D-HST and CAN- 
DELS Treasury surveys. The surface brightness profiles were 
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found to be well represented by Sersic profiles with average 
indices evolving from n ^ 1 (exponential) to n ^ 3 (more cen¬ 
trally concentrated) from z ^ 2.5 to z ^ 0, with ~ 1 dex associ¬ 
ated increase in average stellar masses and a factor of ~ 2 in¬ 
crease in galaxy sizes. These observational trends are close to 
the evolution found in our fiducial galaxy simulation for z < 2, 
which develops a more centrally concentrated mass profiles at 
late times via global disk instabilities, with roughly a Sa-Sab 
z = 0 Hubble type. 

At lower stellar masses (M* < 10'° M.q), the stellar frac¬ 
tion M*/M 2 oo decreases rapidly (see Figure 2), which means 
that these galaxies are more stable against global disk instabil¬ 
ities which prevents formation of central mass concentration. 
Indeed, bulgeless disks are observed to be more abundant at 
these masses (Nair &. Abraham 2010). However, they are also 
more susceptible to disturbances due to injection of momen¬ 
tum and energy by stellar feedback, which can result in irreg¬ 
ular, rather than thin disk, morphology. 

4.2. Comparison with previous studies 

A number of recent studies have emphasized the outmost 
importance of stellar feedback in shaping basic properties 
and morphologies of galaxies. Following initial findings of 
Navarro & White (1994), Navarro & Steinmetz (2000) and 
Abadi et al. (2003) have conhrmed the fact that when gas is 
not removed by winds, or when galaxies are allowed to form 
stars at a high rate during early epochs of evolution, the result¬ 
ing z = 0 galaxies are too compact and have specific angular 
momentum that is too low compared to observations of late 
type galaxies. 

Studies by Okamoto et al. (2005), Governato et al. (2007), 
Zavala et al. (2008) and Scannapieco et al. (2008) have con¬ 
hrmed that the problem arises from the loss of angular mo¬ 
mentum the cold, condensed baryons experience after they are 
accreted in clumps of stars in gas and that the loss of angular 
momentum by baryons is mitigated by increasing the efficacy 
of stellar feedback. 

Efficient feedback prevents baryons from condensing and 
accumulating in dense central “blobs” and delays star forma¬ 
tion until lower redshifts, when the gas is able to settle into 
a gaseous disk without signihcant loss of angular momentum 
(Mailer & Dekel 2002). In particular, Zavala et al. (2008, see 
also Sales et al. 2010, Pedrosa et al. 2014, and Genel et al. 
2015) have pointed out that agreement with observations re¬ 
quires that baryons that end up in late-type galaxies have spe- 
cihc angular momentum comparable to that acquired via grav¬ 
itational torques by their host dark matter halos, as envisioned 
in the classical analytic models of disk formation (Fall & Efs- 
tathiou 1980; Dalcanton et al. 1997; Mo et al. 1998, although 
see Dutton & van den Bosch 2012 who concluded that ob¬ 
served angular momenta of disk galaxies are somewhat lower 
than expected for their DM halos). The same conclusion can 
be reached by comparing sizes of galaxies to the expectations 
of such analytic models (see Kravtsov 2013). 

The importance of efficient feedback in setting galaxy an¬ 
gular momentum, size, and morphology was confirmed by a 
number of recent studies (Sales et al. 2010; Scannapieco et al. 
2012; Brook et al. 2011; Aumer et al. 2013, 2014; Ubler et al. 
2014; Christensen et al. 2014b, a; Genel et al. 2015; Murante 
et al. 2015). Thus, the role of feedback and the importance of 
proper modelling of the star formation-feedback loop are now 
hrmly established. The exact mechanisms by which feedback 
affects galaxy properties are nevertheless still a subject of de¬ 


bate. 

One possible effect of feedback is to self-regulate star for¬ 
mation, making it inefficient and maintaining galaxies gas 
rich. Disks in gas rich galaxies are much more resilient during 
mergers and a thin disk can reform even after major mergers 
(Robertson et al. 2006). This effect is likely significant at 
high redshifts when the age of the universe is comparable or 
smaller than the typical gas consumption time scale in galax¬ 
ies. However, by itself self-regulation cannot result in forma¬ 
tion of realistic disks as the accumulating gas eventually turns 
into stars, and without gas removal via winds, galactic stellar 
masses will end up too high to be compatible with observa¬ 
tions. 

Zavala et al. (2008) and Scannapieco et al. (2008) argued 
that the role of feedback was primarily to prevent gas conden¬ 
sation into dense clumps in low-mass dark matter halos and 
delaying gas accretion onto disk until later epochs, as argued 
by Mailer & Dekel (2002). Brook et al. (2011) and, more 
recently, Ubler et al. (2014) additionally argued that stellar 
feedback can promote disk formation by preferentially eject¬ 
ing low angular momentum gas at early times (z > 1), which 
can accumulate due to disk instabilities or angular momen¬ 
tum loss in mergers. Genel et al. (2015) also emphasized 
the importance of stellar feedback driven winds with large 
loading factors for the formation of angular momentum rich 
disks. They concluded, however, that AGN feedback via ra¬ 
dio bubbles in their implementation leads to a net loss of an¬ 
gular momentum, concluding that AGN feedback may help 
explain low angular momentum of spheroidal stellar systems. 
Overall, Genel et al. (2015) concluded, as did another re¬ 
cent simulation studies by Aumer et al. (2013) and Eiacconi 
et al. (2015), that diverse galaxy morphologies observed in 
their Illustris simulations did depend on the merger history of 
galaxies (as argued since the earliest galaxy formation sim¬ 
ulations Navarro & Benz 1991). Both disk and spheroidal 
galaxies in their simulations had similar specific angular mo¬ 
mentum at high z comparable to that of their dark matter ha¬ 
los. Spheroidal galaxies, however, lost significant fraction of 
their angular momentum in major mergers. Sales et al. (2012), 
however, argued that alignment of angular momentum of ac¬ 
creting gas was a more important factor in setting galaxy mor¬ 
phologies than the merger history. In reality, it is likely that 
both mergers and angular momentum alignment of accreting 
gas affect galaxy morphology and structure (see, e.g., Aumer 
et al. 2014). 

Our conclusions on the high importance of stellar feedback 
in setting properties of galaxies are in broad agreement with 
results of the studies discussed above. Our study have con¬ 
tributed to our continuing quest to understand effects of feed¬ 
back in several ways. Eirst, our simulations implement a num¬ 
ber of state-of-the-art prescriptions for star formation and stel¬ 
lar feedback, such as molecular hydrogen based star forma¬ 
tion and and implementation of both “early” (pre-supernova) 
and supernova momentum and energy injections. Most im¬ 
portantly, our simulations reach a resolution of ^ 50- 100 pc 
in the ISM, allowing for the implementation of star formation 
and feedback prescriptions to operate on these scales, hence 
more closely resolving the internal processes that shape galax¬ 
ies structurally. 

Large-scale galactic winds then emerge (if they do) self- 
consistently due to coherent effects of local feedback from 
multiple star formation sites. This is in contrast to many re¬ 
cent studies, in which large-scale winds are launched explic- 
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itly as part of the overall feedback prescription, often accom¬ 
panied by turning off the interactions between the wind and 
the galaxy’s ISM and circumgalactic gas. For example, at 
z = 0, gravitational forces in the Illustris simulation analyzed 
by Genel et al. (2015) are softened on scales of ^ 710pc, com¬ 
pared to ^ 75 pc in the our current work. Similar difference 
in resolution exists with the EAGLE and the Magneticum 
Pathfinder simulations (Schaye et al. 2015; Teklu et al. 2015). 
The flip side, however, is that the computational expense de¬ 
manded by such high resolution allowed us to simulate only a 
handful of galaxy models. 

Our results indicate that in addition to the role of feedback 
in ejecting low angular momentum gas pointed out in previous 
studies, feedback may also play a role in limiting disk mass, 
thereby suppressing disk instabilities and preventing associ¬ 
ated loss of angular momentum in the first place. 

At the same time, our results show that efficient feedback 
that suppresses galaxy mass to the observed levels or below 
is no guarantee of forming a galaxy with realistic properties. 
Thus, feedback in our ALL_Efb_e001_5ESN run main¬ 
tains a low star formation rate and stellar mass, but the galaxy 
that forms is highly irregular and has size that is too large 
compared to observed galaxies. This is reminiscent of the 
recent results of Roskar et al. (2014), who demonstrated the 
excessively violent energy and momentum input can lead to 
overheating of gaseous disks, even as the overcooling problem 
is alleviated (see also Stinson et al. 2013b). This overheating 
then produces disk galaxies that are much thicker than ob¬ 
served. Although gas stirring by feedback may actually play 
a role in setting the disk thickness (Bird et al. 2013), these re¬ 
sults illustrate that stirring effects should be limited as not too 
thicken the disk excessively. 

Overall, our results show that basic properties of galaxy 
forming in a given halo are highly sensitive to the implemen¬ 
tation and parameters of star formation and stellar feedback, 
consistent with findings of Okamoto et al. (2005) and Scan- 
napieco et al. (2012). This inherent sensitivity to details of 
the star formation-feedback loop indicates that specific re¬ 
sults, for a particular set of prescriptions, should be viewed 
with healthy skepticism, regardless of how well they repro¬ 
duce particular set of observations. 

Our results published in a companion paper (Liang et al. 
2015), in which we compared properties of circumgalactic 
gas in our simulations to observations, also highlight this 
point. Although our fiducial run does result in many realistic 
galaxy features, including stellar mass, morphology, disk-to- 
bulge ratio, rotation curve shape, disk thickness (see sections 
above and results in Agertz & Kravtsov 2015), it fails spectac¬ 
ularly in reproducing observed properties of circumgalactic 
gas around galaxies of similar mass. In particular, its gaseous 
halo is too hot and is almost devoid of any extended warm gas 
component, the existence of which is revealed by ubiquitous 
absorption lines of low-ionization energy ions, such as Mg II. 

The model with boosted energy per supernova 
(ALL_Efb_e001_5ESN ) produces a circumgalactic 
medium closer to observations due to its continuing large- 
scale outflows to low z, fails to produce a realistic galaxy. 
Although the discrepancy with observed CGM properties 
may be due, in part, to low spatial numerical in galactic halo 
regions (Muzahid 2014; Crighton et al. 2015) or missing 
physical processes such as cosmic rays (EnBlin et al. 2007; 
Booth et al. 2013; Salem & Bryan 2014), the fact that a 
comparison with the CGM properties provides “orthogonal” 
constraints illustrates the importance in comparing as wide 


a range of galaxy properties as possible (see also Brook 
et al. 2012b), in particular with the properties of gaseous 
component both within galaxies and in their surrounding 
halos. 

Einally, our results show that changes of parameters related 
to the star formation and stellar feedback prescriptions can re¬ 
sult in variations of galaxy size, ri/ 2 , by more than an order 
of magnitude, while the final specific angular momentum in 
all baryons (Jstar+gas, Eigure 6) varies only by a factor of two. 
This shows that a comparison to both galaxy sizes and an¬ 
gular momenta is warranted, as they provide complementary 
information. Size evolution as a function of galactic mass, in 
particular, will likely to be a useful and stringent constraint 
on galaxy formation models (van Dokkum et al. 2013; Aumer 
et al. 2014). 

5. CONCLUSIONS 

In this study we use cosmological zoom-in simulations of 
galaxy formation in a MW-sized halo to investigate the evo¬ 
lution of galaxy sizes, morphologies and angular momenta in 
runs with different parameters of the star formation - feed¬ 
back loop implementation. In particular, we focus on three 
runs selected from a larger set of seven runs reported in Agertz 
& Kravtsov (2015), which employ the model for star forma¬ 
tion and feedback described in detail in Agertz et al. (2013). 
In these three runs the star formation and feedback prescrip¬ 
tions are kept the same, but the local star formation efficiency 
per free fall-time is varied by a factor of ten, from eff = 0.1 in 
the fiducial run to eff = 0.01 in the two low star formation effi¬ 
ciency runs. One of the Eff = 0.01 runs has otherwise identical 
parameters to the fiducial run (inefficient feedback), while in 
the other the energy released by supernova is increased by a 
factor of five relative to the fiducial value of 10^* ergs. Our 
main results and conclusions can be summarized as follows. 

1. One of the runs - the fiducial run with the star forma¬ 
tion efficiency of eff = 0.1 - generates vigorous large- 
scale feedback-driven outflows during early stages of 
galaxy evolution and reproduces the basic properties of 
late-type galaxies: morphology dominated by a kine¬ 
matically cold disk, size of the stellar distribution, and 
specific angular momentum. In addition, the stellar and 
cold gas surface density profiles in this model is in good 
agreement with profiles of late type galaxies from the 
THINGS sample (Leroy et al. 2008). 

2. The stellar half-mass radius in the fiducial simulation 
evolves approximately along the observed z = 0 ri/ 2 - 
M* relation, and agrees well with the rj /2 - evolu¬ 
tion for MW-like progenitors traced back to z 2.5 by 
van Dokkum et al. (2013) and Patel et al. (2013). It also 
roughly follows the approximately linear ri/ 2 “^ 200 c 
relation (Kravtsov 2013), especially if we correct for 
effects of pseudo-evolution on the halo mass and ra¬ 
dius. The evolution along the observed ri/ 2 -R 200 c re¬ 
lation indicates that the specific angular momentum of 
the galaxy in this simulation is roughly comparable, or 
somewhat smaller, than the angular momentum of the 
dark matter halo throughout galaxy evolution. 

3. Our results indicate that in the fiducial model, 
feedback-driven outflows delay the bulk of star forma¬ 
tion to z ^ 2. After the last major merger, the galaxy en¬ 
ters an epoch of disk formation, with significant growth 
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due to accretion from the hot gaseous halo. Here the 
specific angular momentum of the stars ( 7 *) increases 
by an order of magnitude towards values compatible for 
galaxies of roughly Sa-Sab Hubble type. 

4. Although the implementation of star formation and 
feedback in the run with inefficient local star forma¬ 
tion (Cff = 0 . 01 ) is identical, large-scale outflows are ab¬ 
sent, indicating that energy and momentum injection is 
too weak. The resulting galaxy suffers from the well- 
known overcooling and angular momentum problem, 
and does not resemble a realistic galaxy. We thus con¬ 
firm conclusions of many previous studies that efficient 
feedback is critical to produce realistic late-type galax¬ 
ies. 

5. Our results indicate that one of the roles of stellar feed¬ 
back is to maintain stellar mass-halo mass fractions 
near or below the critical threshold for global disk in¬ 
stabilities. We argue (see Section 4.1) that the failure 
of many previous simulations in predicting the exis¬ 
tence of late-type galaxies, with moderate or low bulge- 
to-disk ratios, in large part is related to the failure of 
maintaining the mass of the gaseous disk below such 
threshold, which leads to strong bars and spiral insta¬ 
bilities that efficiently and quickly channel gas towards 
the central regions. 

6 . At the same time, we show that feedback can be too 

efficient and prevent the formation of a dominant, kine¬ 
matically cold disk component. The run with = 0.01, 
but boosted energy per supernova, produces a diffuse 
irregular galaxy with size much larger than a typical 
galaxy of that stellar mass, but with a rather low angu¬ 
lar momentum. The galaxy in this run deviates from 
the observed relation, which illustrates that 

increasingly abundant observations of the evolution of 
galaxy sizes will be a useful and stringent constraint on 
galaxy formation physics (see also Brooks et al. 2011; 
Aumer et al. 2014). 

7. Although the galaxy produced in the fiducial run is far 
more realistic than the galaxy produced in the run with 
Eff = 0 . 01 , but boosted energy per supernova, our recent 
comparison of properties of circumgalactic gas in these 
simulations to observations presented in Liang et al. 
(2015) shows that the fiducial run fails to reproduce the 
observed CGM properties, while the run with boosted 
Esjsi, leading to an unrealistic galaxy, is a much close 
match to the observations. This illustrates the impor¬ 
tance of testing star formation and feedback prescrip¬ 
tions on the full array of available data. 
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